Non-uniform phases in metals with local moments 
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The two-dimensional Kondo lattice model with both nearest and next-nearest neighbor exchange 
interactions is studied within a mean-field approach and its phase diagram is determined. In par- 
ticular, we allow for lattice translation symmetry breaking. We observe that the usual uniform 
inter-site order parameter is never realized, being unstable towards other more complex types of 
order. When the nearest neighbor exchange Ji is ferromagnetic the flux phase is always the most 
stable state, irrespective of the value of the next-nearest-neighbor interaction J2. For antiferromag- 
netic Ji, however, either a columnar or a flux phase is realized, depending on conduction electron 
filling and the value of J2. 
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The nature of the various magnetic phases of heavy 
fermion compounds has been the focus of attention over 
the years. Most of the analysis is based on the celebrated 
paradigm of Doniach, who conjectured a phase diagram 
consisting of two possible phases, one paramagnetic and 
another exhibiting long range antiferromagnetic order, 1 
The driving mechanism behind this phase diagram is 
the competition between the Kondo effect which favors 
paramagnetism and is dominant at strong exchange cou- 
pling, and the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction al 4 : 5 which dominates at weak coupling and 
can lead to antiferromagnetism. Particularly interesting 
is the quantum phase transition which separates the two 
phases at zero temperature and which can be accessed by 
tuning the exchange interaction between local moments 
and conduction electrons through external or chemical 
pressure. This quantum critical behavior has been in- 
tensively studied experimentallySiS but a complete the- 
oretical description is still lackingiS*ifl*ii^Siiiid 

Despite the appealing simplicity of the Doniach 
phase diagram the possibility of the existence of 
other kinds of phases remains. Among these we 
should mention inhomogeneous magnetic order^i 5 - or- 
bital antiferromagnetismifi and dimerizationiii The last 
possibility has been given strong numerical support in 
the one-dimensional case at quarter conduction electron 
fillingjii It was ultimately ascribed to the long-ranged 
RKKY interaction between localized spins Although 
dimerization is an oft-encountered instability in one di- 
mension, its presence in higher dimensions is less fre- 
quent. There is some (controversial) evidence in favor 
of its existence in the frustrated two-dimensional Heisen- 
berg model with both nearest- and next-nearest-neighbor 
interactions pi2iiS*2Si21 However, the long-ranged nature of 
the RKKY interaction makes its appearance more likely 
in metallic systems with local moments. Motivated by 
this, the aim of the present study is to look for dimeriza- 
tion in particular and other forms of order with broken 
lattice translational symmetry in general in higher di- 
mensional models of heavy fermion materials. 

The co-existence of magnetic inter-site correlations and 



the Kondo effect has been investigated before using mean 
field calculations. Usually, two order parameters are con- 
sidered: one describing the local correlations generated 
by the Kondo effect and the other connected to non-local 
inter-site correlations. If the inter-site correlations break 
spin SU(2) symmetry, there is a competition between 
Kondo singlet formation and magnetic ordering of some 



type- 



Alternatively, the tendency for Kondo com- 



pensation can be analyzed in a scaling approach i 25 i 2 fii 2 T 
On the other hand, if the inter-site correlations do not 
break spin SU(2) symmetry, there may be the formation 
of some kind of spin liquid statei 2 ^' 2 ^^ Fluctuations be- 
yond mean field have also been considered in connection 
with the quantum critical behavior of the system^i*^ 
In this study, we have allowed for the emergence of bro- 
ken lattice translational symmetry in the non-local cor- 
relations, without a broken SU(2) symmetry. We have 
studied the effects of conduction electron filling and both 
nearest- and next-nearest-neighbor exchange interactions 
on the possible phases of the Kondo lattice model in two 
dimensions. The inclusion of further-neighbor interac- 
tions is intended to partially incorporate the long-ranged 
nature of the RKKY interaction between localized mo- 
ments. We have found that the usually assumed uniform 
state is unstable throughout the phase diagram towards 
either columnar or flux phase order £2/21 We have also 
studied the temperature dependence of the order param- 
eters. They do not seem to differ much from the uniform 
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The Kondo lattice Hamiltonian is given by 



H K = 



E 



3,aP 



(1) 



where ek is the band dispersion, fi is the chemical po- 
tential, Cj a and c kcr are conduction electron annihila- 
tion operators in real (Wannier) and reciprocal spaces 
respectively, Sj is a localized spin- 1 ; operator, and <J a j} 
are Pauli matrices. In addition to the above terms we 
also include Heisenberg-like interactions between nearest 
neighbor and next-nearest neighbor localized spins in an 
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attempt to partially capture the long-ranged nature of 
the RKKY interaction. Hence, the full Hamiltonian can 
now be written as H = Hk + Hh, where 



H H = JiJ2 S > 



Sk + J: 



E 



S; • S„ 



(2) 



where (jk) and ((lm)) denote nearest-neighbor and next- 
nearest neighbor sites, respectively. In this work, we con- 
sider both ferromagnetic and antiferromagnetic exchange 
interactions. The spin operators can be expressed in the 
usual Abrikosov pseudo-fermionic representation 
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where a constraint of single /-electron occupancy is im- 
plied. The mean-field Hamiltonian can be written by ex- 
pressing the spin fields in terms of the above /-fermionic 
operators and defining the following three order param- 
eters 



Xi 



4>ja - 2 (djafja + -f}a c ja) > 

ft f i ft f 
JjaJka 1 Jka J ja 

ft f I ft f 
J la J ma 1 J ma •> U 



Xjka 



(3) 
(4) 
(5) 



where j and k are nearest neighbor sites and I and m 
denote next-nearest neighbors. We will focus on SU(2) 
invariant states, hence none of the order parameters will 
depend on a (<f> ja = <f>j, Xjka = Xjk, xLa = Xim)- We 
can write down the mean field Hamiltonian as 



-ko-^kcr 



E fjafja 



H M f = E( £k ~^) c k 

k,<r 

- 2.7* J> (s>4* ' !l - 

j," 

"JiE (xjkflj^+Ti 

{jk),a 

E \X'lmflafm,a + H. C 

((lm)),a 



i^ fe i 2 



{(lm)) 



We will focus on a two-dimensional tight-binding disper- 
sion relation for the conduction band 



D 



(cos k x a + cos k y a) 



(7) 



where D is the half bandwidth and a is the lattice param- 
eter. The chemical potential is determined by the con- 
duction electron density n through jfJ2ka( c t.a c ka) = n 
(N is the number of lattice sites) and E is a Lagrange 
multiplier used to impose the /-electron single occupancy 
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Figure 1: Schematic picture illustrating the various possible 
phases: (a) Uniform phase, (b) Dimer phase, (c) Columnar 
phase, and (d) Flux phase. 



constraint on the average A Y^Uaifkafkcr) ~ The f ree 
energy can be written as 

F = -2T H 1 + e ~ K/T } + ( E o- V n ) N 

k,a=± 

+ 4J A ^|0/ 2 + 2J 1 E| X , fe | 2 + 2J 2 £ \ X ' lm \U 

J (jk) ((lm)) 

where T is the temperature and are the non- 

interacting bands of the mean field Hamiltonian JfjJ • 

In this work, we have chosen energy and length units 
such that both D and o are equal to 1. When trans- 
lational invariance is not broken, <f)j — <fi, \jk — X an d 
Xim = which we will henceforth call the uniform state. 
In addition to the uniform case, we also consider the 
dimerized state with dimers along the x-axis, the colum- 
nar phase and the flux phase i^^ 4 In all cases, and x'i m 
are taken to be uniform. The four phases are described 
as follows (see Fig.^J: 

(a) Uniform: All x' s are real and equal. Lattice transla- 
tion symmetry is not broken. 

(b) Dimers: x ls zer0 f° r bonds along the y-direction 
whereas for bonds along the x-direction we have 



X 

Xjk = 2 



This phase has broken lattice translational and rotational 
symmetries. Another state with the x and y directions 
interchanged is degenerate and equivalent to this one. 
(c) Columnar: x is uniform and equal to x f° r bonds 
along the y-direction whereas for bonds along the re- 
direction we have 



X 

Xjk - 2 



i + (-iy 



This phase also has broken lattice translational and ro- 
tational symmetries. By interchanging x and y direc- 
tions we once again get another degenerate and equiva- 
lent state. 

(d) Flux Phase: All of x' s are equal in magnitude 
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n=0.9, J =0.5 




Figure 2: Phase diagram of the Kondo-Heisenberg model at 
T = as a function of the nearest-neighbor exchange Ji and 
the conduction electron filling n. The next-nearest-neighbor 
coupling J 2 = and the Kondo coupling Jk =0.5. 



Figure 3: Phase diagram of the Kondo-Heisenberg model at 
T = with nearest-neighbor and next-nearest neighbor ex- 
change at n — 0.9 and Jk = 0.5. 



but may have imaginary phases. The specific choice 
of these phases is not gauge- invariant i*-- However, the 
flux through a plaquette is a gauge-invariant quantity. It 
is given by the phase of the oriented plaquette product 
I"[ = X12X23X34X41 "~ We consider the case depicted in 
Fig. [Hd), in which = ±, staggered between adjacent 
plaquettes, corresponding to fluxes of ±7r. This choice 
can be realized by the following gauge choice 

Xjk = \x\ 
for bonds along the y-axis and 

Xjk = (-l) j i\x\ 

for bonds along the x-axis. The bonds now being complex 
have a definite direction which is shown in Fig.QJd). 

We first study the phase diagram of the model at T = 
by varying J\ and n while keeping J 2 = 0. The Kondo 
coupling is kept at Jk = 0.5. This is shown in Fig |2 
We consider both ferromagnetic and antiferromagnetic 
values of J\. The first thing to notice is the instability 
of the uniform state, which is usually assumed, towards 
other forms of order. For antiferromagnetic coupling be- 
tween the local moments the columnar and flux phases 
are the most stable, the latter occurring only for suffi- 
ciently large J\. However, when J\ < (ferromagnetic 
coupling), the flux phase is the most stable, irrespective 
of the filling and the value of J\ . The transition between 
flux and columnar phase is first order. 

We now proceed to investigate the influence of the 
next-nearest-neighbor coupling J2 between the local mo- 
ments, still at T = 0. We studied the phase diagram at 
n = 0.9 and n = 0.4 (Figs. and QJ respectively). We 
have allowed for both antiferromagnetic and ferromag- 
netic couplings between next-nearest neighbors. For fer- 
romagnetic Ji , the flux phase is dominant irrespective of 
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Figure 4: Phase diagram of the Kondo-Heisenberg model at 
T = with nearest-neighbor and next-nearest neighbor ex- 
change at n — 0.4 and Jk = 0.5. 



the value of J2. For antiferromagnetic J\, columnar and 
flux phases share the parameter space. For J\ < 0.24, 
only the columnar phase is realized. For higher values 
of Ji , a flux phase can appear if the conduction electron 
filling is large enough, as shown in Fig. | At n = 0.4, 
on the other hand, the most stable ground state is deter- 
mined solely by the sign of J\ , irrespective of the value 
of J2. In this ferromagnetic Ji favors the flux 

phase, whereas an antiferromagnetic J\ leads to a colum- 
nar phase. Again, the phase boundary between flux and 
columnar phases is a first order line. 

In Fig. El we show the filling dependence of the order 
parameters <f> and x at T = 0, for Ji = 0.2, J 2 = 0, and 
Jk = 0.5. In this case, the system is always in a columnar 
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Figure 5: Filling dependence of the order parameters 4> and 
X at T = in the Kondo-Heisenberg model with nearest- 
neighbor exchange only (Ji = 0.2) and Jk = 0.5. 

phase (see Fig.EJ. There is a clear competition between 
the two types of order, the Kondo effect (<j>) becoming 
more predominant as the system approaches half-filling. 
Of course, this competition is analogous to the one pre- 
dicted by Doniach between a tendency to form to local 
singlets (0) and another one to lock localized spins into 
some kind of order. Our mean field Ansatz is able to cap- 
ture this competition. The predominance of the Kondo 
effect as the system approaches half-filling is due to an 
enhanced density of states in that region providing more 
conduction electron states to quench the local moments. 
By contrast, note that, for the same parameters of Fig.0 
the uniform order parameter has a much more reduced 
value and does not compete with the Kondo effect at the 
mean field level (see Fig. 3 of Refill) . 

In addition, we have also studied the temperature de- 
pendence of the order parameters for J\ = 0.2, J2 = 0, 
and Jk — 0.5. The temperature dependence has been 
plotted for n = 0.8 in Fig. Eland n = 0.4 in Fig. □ In 
both cases, the columnar phase is the most stable from 
T = up to the transition temperature. Although the 
two dependences are different, both order parameters <f> 
and x disappear at the same critical temperature. This 
same simultaneous disappearance of order had been ob- 
served in previous studies of the uniform phase for similar 
values of the exchange couplings^ Although the finite 
temperature phase transition triggered by x could be re- 
alized in real systems, the vanishing of <j) is an artifact of 
the mean field treatment^ 

Let us now pause to compare our results with previous 
studies. A mean field Ansatz of the form considered here 
has been investigated before ; 2 fi i 2 Mft without allowance for 
broken lattice translation symmetry. An important con- 
clusion of our results is that the uniform state consid- 
ered in these references is never stable. Refs. 
on the other hand, do consider the effects of both near- 
est and next-nearest neighbor couplings between local- 
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Figure 6: Temperature dependence of the order parameters <j> 
and x (columnar phase) at n = 0.8 in the Kondo-Heisenberg 
model with nearest-neighbor exchange only (Jj = 0.2) and 
Jk = 0.5. 
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Figure 7: Temperature dependence of the order parameters <j> 
and x (columnar phase) at n — 0.4 in the Kondo-Heisenberg 
model with nearest-neighbor exchange only (Ji = 0.2) and 
Jk = 0.5. 



ized spins. Their treatment of the Kondo effect, how- 
ever, is confined to a scaling analysis, which breaks down 
below the Kondo scale. Our self-consistent treatment 
of the Kondo effect, by contrast, is able to reach deep 
into the Kondo singlet formation regime and thus offers 
a better treatment below the Kondo scale. Finally, no 
comparison has been attempted with the mean field free 
energies of phases with conventional long-range magnetic 
orde n 22 i 2 i?i 24 This would determine the region of stabil- 
ity of these non-uniform phases. We leave this for future 
studies. 

In conclusion we have studied the mean field phase di- 
agram of the two-dimensional Kondo-Heisenberg model 
with both nearest- and next-nearest-neighbor exchange 
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interactions for various values of doping, temperature 
and coupling constants. We have observed that the uni- 
form state solution is unstable towards lattice transla- 
tional symmetry breaking for any value of the exchange 
constants. Depending on the values of J\ , J2 and filling n, 
the system realizes either a columnar or a flux phase. The 
flux phase is always stabilized by a nearest-neighbor fer- 
romagnetic exchange between localized spins. When this 
coupling constant changes sign, however, both columnar 



and flux phases can occur, the latter being favored at 
large J\ and n and the former appearing at small J\ and 
low fillings. 
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